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ABSTRACT: The Dirac equation is solved in the Einstein- Yang-Mills background found 
by Bartnik and McKinnon. We find a normalizable zero-energy fermion mode in the s-wave 
sector. As shown recently, their solution corresponds to a gravitational sphaleron which 
mediates transitions between topologically distinct vacua. Since the Bartnik-McKinnon 
solution is unstable, it will either collapse to form a black hole or radiate away its energy. 
In either the Chern-Simons number of the configuration changes, there will be an 

accompanying anomalous change in fermion number. 
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1. Introduction 

Neither tiie Yang-Mills nor Einstein field equations admit static finite-energy, non- 
singular solutions. However, such particle-like solutions have been found by Bartnik and 
McKinnon in the combined Einstein- Yang-Mills theory. These solutions were later 
interpreted as sphalerons, that is, static saddlepoint solutions lying at the top of an energy 
barrier in field configuration space separating vacua with different Chern-Simons number 
0. This, in particular, accounts for their instability. If the sphaleron is perturbed, it will 
either radiate its energy to infinity or collapse to form a black hole. Since either process 
involves a change in Chern-Simons number, one expects an equal anomalous change in 
chiral fermion number. In this paper we initiate a study of this problem. In Section 2, the 
Bartnik-McKinnon solution and its interpretation as a sphaleron is reviewed. In Sections 
3-5, the Dirac equation in the background of the static Bartnik-McKinnon sphaleron is 
analyzed. We find a zero energy bound state in the s-wave sector. In Section 6, a Higgs 
field is included. In Section 7, the conformal invariance of the Dirac equation is exploited 
to prove a general no-hair theorem for fermions. 



2. Bartnik-McKinnon Solutions 

In |jl|] a class of spherically symmetric particle-like solutions to the Einstein- Yang-Mills 
equations with SU{2) gauge group were found numerically. It was shown that there are 
an infinite sequence of solutions (A^(r), i?^(r), iir^(r)) for the metric and gauge field 

1 — K 

Ai = — e,jfcn^r'= * (2.2) 

2gr 

where g is the gauge coupling constant, r* are the generators of SU{2), and 
fi = (sin 9 cos ^, sin 9 sin (p,cos 9). has n nodes and satisfies the boundary con- 
dition Kn{0) = 1 and Kn{oc) = (—1)". At large distances, the metric approaches 
Schwarschild and the gauge field strength decays as with Ai approaching the pure 
gauge Ai = —^diUU~^ where U = 1 for n even and U = —in ■ r for n odd. Near a 
node r = ro (-ftr(ro) = 0), the gauge field corresponds to a, P = 1/g Dirac monopole, 
and the metric to extremal Reissner-Nordstrom. The masses of the Bartnik-McKinnon 

* This form of the gauge field is gauge equivalent to that used in [jl]. 



1 



solutions are proportional to the only mass scale in the problem, g~^G~^/'^ , where G is the 
gravitational constant. The masses increase with n approaching the mass of the extremal 
Reissner-Nordstrom black hole with magnetic charge, P = 1/g a.s n oo. Heuristically, 
in this limit, fluctuates rapidly approaching its mean zero value corresponding to the 
Dirac monopole. The existence of the Bartnik-McKinnon solutions was subsequently es- 
tablished rigourously [Q, and generalized to include horizons 0. 

Soon after their discovery, it was shown that both particle-like and black hole solutions 
are unstable In fact, the Bartnik-McKinnon solutions correspond to a gravitational 
analog of the sphaleron 00. We recall that sphalerons are static saddle-point solu- 
tions that lie at the top of an energy barrier separating topologically distinct vacua. For 
asymptotically flat spacetimes, the ADM energy provides a positive deflnite mass func- 
tional on fleld conflguration space. The zero energy "vacua" are flat spacetime metrics 
and pure gauge Yang-Mills conflgurations Ai = —^diUU~^. Each pure gauge yields a 
map from space into the group manifold of SU{2). Demanding U 1, this becomes a 
map U : S^. Since maps with different winding, or Chern-Simons, numbers cannot 

be continuously deformed into one another, paths in conflguration space connecting these 
vacua must enter no n- vacuum regions, and therefore by the positive deflniteness of the mass 
functional must traverse an energy barrier. * Each such path has a maximal energy conflg- 
uration. Among these conflgurations the one with minimum energy will be an extremum 
of the mass functional and therefore correspond to an unstable saddlepoint solution. By 
considering paths connecting the vacuum sector with Chern-Simons number zero and the 
sector with Chern-Simons number one, one obtains the Bartnik-McKinnon solution Ki. 
One can see by symmetry that its Chern-Simons number is 1/2. Repeating this procedure, 
with paths connecting Ki and its gauge equivalent solution with Chern-Simons number 
— 1/2, the K2 solution with zero Chern-Simons number is obtained. Continuing in this 
way, the entire Bartnik-McKinnon sequence is produced. This procedure also yields the 
black hole solutions if one considers asymptotically flat metrics with flxed horizon area. 

As the flelds evolve dynamically along a path in conflguration space connecting 
vacua with different Chern-Simons numbers and passing through the sphaleron, anoma- 
lous fermion production will occur The total number of chiral fermions created equals 

* We shall regard vacua which are related by so called large gauge transformations as distinct. 
One can, if one wishes, identify them, in which case the paths we refer to below become closed 
non-contractible loops in this identified configuration space. 
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the change in the Chern-Simons number of the field configuration. As one passes between 
topologicaUy distinct vacua, a fermion state will emerge from the negative energy Dirac 
sea, enter the discrete spectrum, and finally merge with the positive energy continuum. 
By symmetry, therefore, there should be a zero energy bound state at the midpoint of this 
path where the sphaleron configuration is located. If perturbed, the sphaleron will fall 
in one of two directions. The Yang-Mills field will dominate, in which case the sphaleron 
will radiate away its energy to infinity. The fermion zero mode would then also escape 
to infinity perhaps after being refiected at the origin. Alternatively, gravity may induce 
the object to collapse to a black hole swallowing the Yang-Mills field, and presumably the 
fermion mode as well. 



3. Dirac Equation 

We now proceed to test this picture by solving the Dirac equation in the sphaleron 
background. As expected from the above discussion, we find a zero-energy bound state 
in the s-wave sector. We also derive the full set of time dependent radial equations. To 
determine more precisely the fate of the zero mode if the sphaleron collapses to a black 
hole would require numerically solving the time dependent equations. Such an analysis 
was done for the Yang-Mills-Higgs sphaleron in . 

The massless Dirac equation in an Einstein- Yang-Mills background for an isodoublet 
fermion is given by 

h''iV^-igA^)<i> = (3.1) 

where is the covariant derivative = — ju'^^'jalb- IJ- and a are tangent and 
spacetime indices respectively and are related by = e", a basis of orthonormal one- 
forms. = w"^ are the associated connection one-forms obeying de" -|- A = 0, and 
7" are Dirac matrices satisfying {7^,7''} = —2r]"-^ with 77^*^ = —1. The metric and gauge 
field are given in (|2.1|) and ( |2.2|) where are now the Pauli spin matrices. In terms of the 



orthonormal basis 

e^ = Adt e'' = Bdr = rdO e"^ = rsinedcj) (3.2) 
the connection one-forms are 

u;^"^ = B-^A'dt uj^'' = B-^de w'^'' = B-isinM(/) u;'^^ = cosM(/) (3.3) 
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where / = ^ . In the chiral representation, the gamma matrices are 



,0 _ / M .,i _ / ctM . / -1 



where a* are Pauh spin matrices. \E' can be decomposed into its left and right chiral 
components \1/ = ^^^^ where i/^LiR) carry two-component Lorentz and isospin indices. 
Since the Dirac equation is massless, the two chiralities decouple. In the following, we 
consider just the right-handed component "0 = i/'Ji- The Dirac equation then becomes 

+ — — (rAV2^) + ^ - igAa ■Ai; = (3.5) 

ot r B or r 



where IpT is the Dirac operator on the unit two sphere and a-h = a* 6*. Eqn. (2^) implies 

K-l 
2gr 



a ■ A = ''^ {n ■ a X r). The conserved inner product is given by 



< ^i|^2 >= j 11^1^2 Br^drdfl. (3.6) 

We now discuss the effect of charge conjugation, C, and parity, P on the sphaleron. 
The action for a Dirac fermion coupled to a Yang-Mills field with arbitrary gauge group 
is invariant under C and P separately. Under C, the fields transform as 



A A'^ — —A* 



(3.7) 



and under P as 



^'(x^^) ^ = 7°*(-x^^) 

A,{x\t) Af = -A,{-x\t) (3.8) 
Ao{x\t) ^ A^ = Aoi-x\t). 

The metric being real and spherically symmetric is invariant under C and P. Because 
C and P interchange chirality, the action for chiral fermions is only invariant under the 
combined CP transformation. Using ( |3.7|) , (|3.8|) , and — r'* = t'^t'^t'^, the gauge field ( |2T^ ) 
transforms as Ai Af^ = t'^Ait'^. The CP invariance of the theory then implies that 
given a solution \1/ to the Dirac equation in the sphaleron background with energy E, 
-j-^^CP jg ^ solution in the same background, but with energy —E. 



4. S'-wave Sector and Zero Mode 



We now solve the Dirac equation ( p.5|) by separating variables. Since the total angular 
momentum K = L + S + T commutes with the Hamiltonian, '0 can be expanded in 
eigenstates of and with eigenvalues k and m. L, S = a/2, and T = f/2 are 
the orbital angular momentum, spin, and isospin. The s-wave sector corresponding to 
/c = m = is spanned by the two states xi ci-nd X2 = ^ ■ ^Xi where 



1 

Xi 



V2 



Jin.- 



(4.1) 



is the hedgehog spinor satisfying {a + t)xi = 0. The action of the various operators 
in the Dirac equation on the two states can be determined from the hedgehog property 
and the spin commmutation relations. The transverse Dirac operator can be written as 
IpT = {2S-L + l){a-n) = {J^ -L^ - + l){a-n) where J = L + S. Since [J^,a-n] = and 
a-fi changes the value of the orbital angular momentum by one, we find that Ipx Xi = — X2 
and IpT X2 = Xi- The operator n ■ a x t appearing in the gauge field term in the Dirac 
equation also interchanges xi ci-nd X2- {n ■ a x r)xi = —2ix2 and {n ■ a x f)x2 = 2ixi- 



Thus, (|3.5D becomes 

^ + ^ + = 

dt dr* r / . ^\ 

dg df AK , ^ ' 

— H / = 

dt dr* r 

where ip = r~^A~^^'^fxi + f~^A~^/'^gx2 and r* is the "tortoise" coordinate satisfying 
^ = A-^B. 

dr 



Eqn. (|4.2|) admits a zero energy bound state of the form 



K 

/ = exp / B— dr , g = (4.3) 



for K = Kn with n odd. For n even, K 1 asymptotically, and therefore / diverges. 
The zero mode ( ^1.3|) vanishes at r = and falls off as 1/r asymptotically. The charge 
density of the wave function from (|3.6|) is given by 47iA-^B\f\^. The density of the zero 
mode ( |4.3|) is peaked in the monopole region, K = 0, which is effectively acting as a 
potential well. For a spacetime with horizon at r = r//, the metric functions behave as A ~ 
(r — rj/)^/^(l + 0(r — r//)) and B ~ (r — rj/)~^/^(H-0(r — r//)). From ( ^.3| ) we observe that 
the wavef unction diverges as (r — rj/)~^/^, and its norm logarithmically. This is perhaps to 
be anticipated in light of the no-hair theorem for fermions. The effective potential within 
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the s-wave sector is found by rewriting the two coupled first order equations as decoupled 
second order equations 



J-aJ,^y(r)f = 0, 9 = -^,+Hf (4.4) 



where 



V{r) = ^+H', H^A-. (4.5) 
or* r 



Near the horizon, the effective potential vanishes as F ~ (r — ri^)^/^. 



5. A; > Modes 

The eigenspaces with eigenvalues k and m are four-dimensional for /c > 1. Since 
and commute, they can be simultaneously diagonalized. Within the four dimensional 
eigenspace, there are two states denoted \k,m, ± > with eigenvalues I — k ±1, and an 
orthogonal two-dimensional space with eigenvalue I = k. In addition, either = (L + S)^ 
or R'^ = {S + T)'^ can be diagonalized, but not simultaneously since they do not commute. 
The eigenstates of are |A;,m, ± > with eigenvalues j — k ± 1/2 and \k,m,j± >, lying 
in the I = k subspace, with eigenvalues j = k ±1/2. The eigenstates of are \k,m,± > 
both with eigenvalues r — 1 and |/c,m,0 > and \k,m,l >, lying in the / = k subspace, 
with eigenvalues r = and 1. The and R^ eigenstates in the I = k subspace are related 
by a rotation matrix of angle ^ with tan^ = sj^^ '■ 

\k, m,j-\. > = {Vk + l\k, m, > +Vk\k, m, 1 >) 

(5.1) 

\k, m, i- > = (—vklk, m, > +\/k + Ilk, m, 1 >). 

V2k + 1 I ' ' / 



The eigenstates when expressed in terms of spherical harmonics take the form: 



Ifc, m, + >= 



1 



^{2k + 2){2k + 3) 



^/{k + m + l){k + m + 2)Y^+^\1, -1 > 



'2{k + m + l){k -m + l)n"^i|l, > +^/{k - m + l){k - m + 2)Y^-^\1,1 > 

1 



>= 



v/2/c(2/c- 1) 



v/(A; - m)(A; - m - -1 > 



+ ^/2{k + m){k - m)Y^^\l, > +^/{k + m - 1) {k + m)Y^_-^ > 



|/c,m,0 >=yfc"'|0,0 > 



|fc,m,l>= , ^/{k -m){k + m + 1)^7^+^11, -1 > 

+mv^yfc"^|l,0 > -v/(A; + m)(A;-m+l)y;"-i|l,l > 



(5.2) 



where |1, 1 >, |1, >, |1, — 1 > are the spin-one triplet of R^. 

We now proceed to express the various operators appearing in the Dirac equation 
as matrices in the basis (|/c,m, + >, |/c,m, j+ >, |/c,m,j_ >, |/c, m, — > ) in which K^^ 

L'^, and are diagonalized. Since a-n commutes with and squares to unity, after 
appropriate normalization of states, one finds that a ■ n\k,m,± >= \k, m,j± >. Using this 
result, one can show that the transverse Dirac operator 1/)t = [J"^ — — S"^ l)a ■ n 
takes the form 

/ -/c-1 \ 
77) I /s + l 

V k / 

Finally, one can determine the matrix form of the operator n ■ a x t appearing in the 



gauge field term in the Dirac equation (3^). Since it is antisymmetric in a and f, it must 
interchange the r = 0, 1 states, and because of the factor n, it changes the orbital angular 
momentum by one. Further calculation shows 



2i 



n ■ a X T 



2k + 1 



( + 1 -^k{k^\) 

-k-\ ^k{k^\) 

y^k(kTT) -k 

V -y/k{k + l) k / 



(5.4) 



Substituting (|5.3|) and ( p.4|) in ( p.5| ), the Dirac equation reduces to four coupled linear 
first-order equations. It appears that higher mode zero energy bound states are forbidden 
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since the wave function diverges at r = 0. Near the horizon, r = rn, A = a{r — rnY^'^i^ + 
0{r — th)) and B = f3{r — rH)~^^'^{'^ + 0{r — th)) while K is constant. Therefore, to 
leading order in r — r^j, the Dirac equation becomes 



(5.5) 



with solutions ip ^ e ^^^(r — ry^)^^'^^^/°' ^1^^ which as before diverge at the horizon. 



6. Higgs Field 

In this section, we consider the effect of a Higgs field. The Dirac equation then becomes 



(6.1) 



where $ = ^V* is the Higgs field in the adjoint representation. We have absorbed the 
Yukawa coupling constant in $. Consider the following spherically symmetric ansatz for 
$ 

0^ = F{r)n\ (6.2) 

Since the two chiralities of the fermion, '^i, and V'fl, no longer decouple, the Dirac equation 
now reduces to four coupled first order equations in the s-wave sector and to eight in the 
higher mode sector. Using the fact that r • nxi = — X2 and r • = — Xi) we find that 
the s-wave equations become 



dfR{L) 
dt 

dgRjL) 
dt 



± 



± 



dgR{L) ^ AK 
dr* r ^^^^^ 

dfR(L) _ AK 
dr* r 



fR{L) 



- lAFgL^R) = 

- iAFh^ji-) = 



(6.3) 



where r(l) = r ^A ^^^/i?(L)Xi + r ^A ^^'^gR(L)X2- There is in fact still a zero-energy 
bound state solution to these equations: 



r K 

fR = -ih = exp / B{ F) dr, gji = gL = 0. 



(6.4) 



Moreover, since the Higgs field causes the wave- function to decay exponentially at infinity, 

is a bound state for both odd and even n. As before, for a space- 
time with horizon the zero mode diverges there. The s-wave equations can be written as 
decoupled second order equations 



R 



Qj-2 Qj.*2 



+ V{r)fR = 0, /l = iJr, qr = iql = 



dfR 
dr* 



+ Hf. 



R 



(6.5) 



where 



V{r) 



dH 

dr" 



H 



r 



(6.6) 



The equations for the higher modes can be found as weU. The matrix associated with the 
operator f ■ n appearing in the Higgs field term is obtained from a ■ nhy exchanging a and 
r and using the fact that the r = 0, 1 states are antisymmetric and symmetric in a and r 
respectively: 



r • n 



( 

—cos 2^ 
sin 2^ 
\ 



-cos 2^ 



sin 2^ 



sin 2^ 



cos 2^ 





sin 2^ 
cos 2^ 




(6.7) 



where tan^ 



7. Fermion No-Hair Theorem 

Assuming that the endpoint of gravitational collapse is a stationary black hole, then 
the uniqueness theorems imply it must be Kerr-Newman. In general, the no-hair theo- 
rems assert that the only external fields that a black hole can generate are those yielding 
conserved charges in the form of a surface integral at infinity. If one attempts to find a 
neighbouring solution with hair, the solution for the perturbation will necessarily diverge 
at the horizon of the black hole (assuming it vanishes at spatial infinity.) We should point 
out that a Yang-Mills black hole such as the one discussed earlier does not unambiguously 
violate the no-hair theorems since the non-linear nature of the field allows one to view it 
and not the black hole as being the source of the external Yang-Mills field. * After all, one 
would certainly ascribe the complicated gravitational field of an accretion disk surrounding 
a black hole to the disk and not to the black hole. The higher multipoles moments of the 
gravitational field can be unambiguously identified with the higher moments of the matter 
density distribution of the disk. The Yang-Mills field, however, does not have well defined 
multipole moments. In particular, attempts to construct a total charge by integrating the 
Yang-Mills magnetic field does not yield a gauge invariant quantity since there is a free 
uncontracted group index. 

The fact that the zero mode ( [4. 3] ) (( |6.4| )) diverges on the horizon is expected from 
these no-hair theorems. In this section, we prove a general no-hair theorem for fermions. 



* In fact, Einstein-Skyrme black holes with Yang-Mills hair have been constructed recently[ll| 
and are stable under linear perturbations [Q. 
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This has aheady been done for Schwarschild in |13] and \T^] By employing the conformal 
invariance of the Dirac equation, we show that aU static solutions in a spherically symmetric 
spacetime diverge at the horizon. Any static, spherically symmetric spacetime may be 
written in isotropic coordinates 

ds^ = -V^dt^ + W^d:>c- (7.1) 

where V and W are functions of |x|. (This ansatz is, in fact, more general than being 
spherically symmetric applying to multi-black hole metrics as well.) Consider the gauged 
Dirac equation 

il^iV^,-igA^)^ -m^ = (7.2) 

where m might depend on a Higgs field and hence on position. In d spacetime di- 
mensions we have the following result: if (\1/, (7^,^, A^, m) is a solution of ( [7. 21 ), then 
(n^~\I/, Q.~'^g^jj, A^, Qm) is also a solution. For the metric ( |7.1D with Aq = 0, this implies 
that solves ( |7.2D with mass term Vm in the ultrastatic optical metric 

ds^ = -dr + -r-T^rfx ■ rfx. (7.3) 

As is the case for no-hair theorems, we are interested in time-independent solutions to the 
Dirac equation. Applying conformal invariance again, but now to the spatial part of the 
metric implies that x = V^^^W^ solves the flat three-dimensional gauged Dirac equation 
with mass term mW. 

For m = Ai = the flat space Dirac equation may be solved in a variety of ways, but 
perhaps the most illuminating from the present conformal viewpoint is to note that in flat 
Euclidean n-space, 

x''=T^,^,...nX^'x^'...x'^ (7.4) 
clearly solves the flat space Dirac equation with m = Ai = provided satisfies 

l'l^T^j...k = 0- (7.5) 

Solutions which vanish at spatial infinity can be obtained from ( [7.4| ) by an inversion. To 
each static solution x{x'^) to the flat space Dirac equation in n spatial dimensions, there is 
an inverted solution ^^x(xV?''^). Applying this to ( |7.4| ), we obtain the general multipole 
solution 

X = ^T,^^,...^y'n^'■■■n'^ (7.6) 
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where n* = jr. Thus, if m = = 0, the general solution in the spacetime which 
decays at infinity is given by 

It is now clear that ^ will blow up at a non-extreme horizon for which V = ^ for r > 0. 
In the extreme case for a regular horizon where W —>■ and V r, we see that ^' still 
blows up for any permitted value of /. 
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